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An  analytical  and  experimental  study  is  being  conducted  to  promote  an  under- 
standing of  the  processes  governing  the  emission  characteristics  of  continuous 
combustion  paver  sources  and  thus  provide  a basis  for  reducing  adverse 
environmental  effects  and  for  controlling  plume  signatures  resulting  from 
aircraft  operations.  The  nodding  dcvclopoent  can,  in  addition,  be  readily 
adapted  to  duns  combustor  and  gas  turbine  combustion.  The  configuration  chosen 
tor  the  study  is  an  opaosed-jet  laboratory  combustor  (OJC).  Three  nunericnl 
programs  arc  used  for  modelPlnc  purposes  ■ The  current  year  has  emphasised  the 
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utility  and  range  of  applicability  of  the  numerical  methods  for  the  case  of 
isothermal  flow.  Eddy  viscosity  madcls  and  boundary  condition  specification 
have  received  the  drearest  emphasis  in  a series  of  systematic  test  of  prediction 
against  experiment.  The  numeriral  methods  are  currently  being  expanded  to 
include  propane  os  veil  ns  methane  oxidation  kinetics.  Evaluation  of  the  numeric 
codes  for  the  ease  of  hot,  reacting  flow  will  be  the  emphasis  of  the  continuation 
year.  Preheat  capability  has  been  added  to  the  experimental  dimension  of  the 
study.  In  addition,  extensive  evaluation  of  chemical  transformation  of  nitrogen 
oxides  while  sampling  combustion  products  is  in  progress. 
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This  second  Interim  Scientific  Report  covers  the 
twelve  (12)  month  reporting  period,  1 May  1976  to  29 
February  1976.  The  first  ten  months  of  the  research 
program  were  structured  to  accommodate  (1)  modification 
and  testing  of  the  experimental  system,  (2)  modification 
of  the  numerical  methods  to  the  geometry  of  the  opposed 
reacting  combustor,  (3)  expansion  of  the  numerical 
methods  to  include  multicomponent  reacting  flow  and, 

(4)  initiation  of  the  cold  flow  numerical  and  experi~ 
mental  studies.  The  current  reporting  period  covers  the 
subsequent  twelve  (12)  month  period  and  includes  (1) 
the  test  of  the  numerical  methods  against  experiment 
for  the  case  of  isothermal  flow,  (2)  expansion  of  the 
numerical  methods  to  Include  propane  oxidation  kinetics 
in  addition  to  methane  oxidation  kinetics,  (3)  expansion 
of  the  experimental  facility  to  Include  air  preheat  capa- 
bility, (4)  generation  of  a broad,  experimental  data  base 
for  the  cases  of  isothermal  and  hot  flow,  and  (5)  an 
in  depth  exploration  of  chemical  transformations  that  may 
occur  in  probes  and  lines  used  to  sample  and  transport 
nitrogen  oxides. 

The  progress  of  the  second  reporting  period  is  present- 
ed and  accompanied  by  a description  of  the  progress 
anticipated  for  the  continuation  year. 


1.  Ih’TkODlJCTIOH 


An  important  category  of  conbuscion  is  designated 
continuous  combustion  because  of  a characteristic  feature 
of  uninterrupted,  chenically  reacting  backmixed  flow.  The 
backmixing  serves  as  a ncchanisn  for  flano  stabilization. 
Partially  or  completely  reacted  species  are  backmixed  with 
incoming  fresh  reactants  to  provide  the  energy  necessary  to 
sustain  the  reaction.  Power  systems  that  operate  on  the 
principle  of  continuous  coabustion  include  gas  turMae 
engines  and  dump  combustors.  As  slown  by  • Se.vtic 

representation  of  a gas  turbine  coabustor  in  Pigu/e  1, 
backmixing  of  inactive  species  can  be  generated  down  stream 
of  the  expao*<on  s-cp  (zone  3)  and  at  the  axis  due  to  swirl 
imparted  to  entering  flow  (zone  A)1. 

1 ^ The  Need  for  Predictive  Methods 

Engineering  deoi^n  of  continuous  combustion  devices  has 
Si  edit  tonally  focuse*’  o:,  v«*rall  thermodynamic  efficiency, 
heat  transfer  r acd  *>owcr  production.  Present  operating 
constraints  (e.g.  exhaust  composition  and  fuel  economy)  have 
prompted  a more  detailed  study  of  the  processes  occurring  within 
the  combustion  zone.  The  more  detailed  studies  can  be  greatly 
facilitated  by  the  application  of  predictive  methods  to  evaluate 
combustion  behavior,  and  the  refined  insight  provided  by  point- 
by-point  predictions  of  concentration,  temperature  nod  velocity 
can  be  especially  valuable.  The  development  of  combustion 
modeling  capabilities  is  also  economically  advantageous  in 
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that  accurate  predictions  of  combustion  performance  would 
reduce  the  need  for  costly  experimental  development  programs. 

The  General  Problem  in  Modeling  Turbulent  flows.  The 
combustion  processes  occurring  in  continuous  systems  represent 
a complex  interaction  of  chemical  kinetics*  heat  transfer, 
mass  diffusion  and  fluid  mechanics.  The  detailed  analyses  of 
combustion-chamber  flows  generally  incorporate  mathematical 
cub-models  of  the  chemical  and  turbulent  transport  processes 
in  a suitable  numerical  framework  for  solving  the  governing 
conservation  equations.  Stemming  from  the  lack  of  a detailed 
knowledge  to  turbulence-combustion  interaction,  a primary 
problem  in  modeling  any  practical  combustion  system  is  to 
adequately  account  for  the  influence  of  turbulent  fluid  motions 
on  the  heat  and  mass  transfer  and  on  the  chemical  reaction  rates. 

The  Additional  Problem  in  Continuous  Combustion.  A 
characteristic  feature  of  practical  continuous  combustion  systems 
is  a zone  of  recirculating  flow  or  backnlxing.  The  zone  serves 
as  a mechanism  for  flame  stabilization.  Partially  or  completely 
reacted  species  are  backmixed  with  incoming  fresh  reactants, 
providing  the  energy  necessary  to  sustain  the  reaction.  As 
a result,  the  zone  of  recirculation  is  a major  factor  in  deter- 
mining the  overall  features  of  the  combustion  flovflcld  and 
must  be  carefully  considered  in  the  application  of  modeling 
methods.  To  account  properly  for  the  zone  of  recirculation 
requires  solution  of  the  appropriate  elliptic  partial  differ- 
ential equations  rather  than  the  less  cumbersome,  more  exten- 
sively investigated  parabolic  equations  applicable  when  addressing 
only  boundary  layer  type  flows. 
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1 . 2 The  Heod  for  F.xpcr incn ta  1 Data 

Experimental  data  ar**  required  to  (1)  evaluate  and 
develop  the  utility  and  range  of  applicability  of  the  pre- 
dictive methods,  and  (2)  to  explore  pollutant  and  plume  for- 
mation in  continuous  combustion  both  prior  and  subsequent  to 
the  availability  of  a suitable  predictive  method. 

The  development  and  evaluation  of  predictive  methods 
depends  upon  the  design  and  conduct  of  a systematic  experimental 
program  to  provide  the  data  necessary  to  assess,  the  numerics, 
boundary  conditions,  initial  conditions,  and  the  submodels 
of  energy  and  mass  transport,  and  chemical  reaction.  The  most 
productive  approach  to  resolving  the  validity,  utility,  and 
range  of  applicability  of  predictive  methods  is  to  perform 
iclevant,  well-controlled  experiments  designed  to  isolate  a 
single  dominant  phenomena,  and  to  perform  in  conjunction  with 
the  experiment,  simultaneous  numerical  simulations.  Progress 
in  continuous  combustion  modeling  is  currently  restricted  by 
the  scarcity  of  detailed  flovfield  data  for  distributed  reac- 
tions with  recirculation.  Measured  profiles  of  chemical  comp- 
osition, temperature,  velocity,  and  turbulence  intensity  in 
backnlxed,  distributed  reactions  are  necessary  to  effectively 
develop  and  evaluate  predictive  methods  for  continuous  combustion. 

In  the  absence  of  a fully  evaluated  and  developed 
predictive  method,  the  experiment  must  guide  the  exploration  of 
pollutant  and  plume  formation  in  continuous  combustion. 

Though  costly  and  cumbersome  in  comparison  to  numerical 
simulation,  experimental  studies  must  currently  be 
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relief  upon  to  provide  insight  Into  the  performance  of 
these  complex  systems.  Nevertheless,  the  parallel  develop- 
ment and  evaluation  of  predictive  methods  can  greatly  enhance 
the  information  derived  from  the  experiment. 

1.?  the  UCI  Program 

Object  1 yes . The  objectives  of  the  UCJ  program  are  directed 
to  clarifying  the  relative  influence  of  those  mechanises  respon- 
sible for  pollutant  production  in  continuous  combustion. 

I 

Specifically  the  objectives  of  the  study  ares 
Near  Term 

o To  investigate  high  temperature  chemical  reactions, 
studying  the  limits  of  combustion  under  those  con- 
ditions typical  of  recirculating,  turbulent  reacting 
flowficlds.  Such  information  is  pertinent  to  com- 
bustion stability  and  use  of  alternative  fuels  in 
turbine,  dump,  and  ramjet  combustors. 

o To  develop  an  understanding  of  pollutant  formation 
in  continuous  combustion  stabilized  by  recirculation. 

Such  information  is  pertinent  to  reducing  environ- 
mental impact  and  controlling  plume  signature. 

o To  develop  and  verify  numerical  methods  and 
associated  submodels  of  turbulence  and  kinetics 
os  applied  to  recirculating  turbulent  reacting 
flowficlds  characteristic  of  turbine  combustion, 

via  a judicious  coupling  of  numerical  methods  to  « 

experiment.  The  goal  is  to  establish  a method  ] 

that  can  be  readily  adopted  to  the  flow  geometry  j 

of  gas  turbine,  dump,  and  ramjet  combustors.  j 

I 

Future  j 

o To  investigate  emission  and  performance  character-  i 

istics  of  alternative  fuels  with  specific  emphasis  ; 

on  particulate  formation  and  flame  luminosi ty . J 

5 
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Approach.  The  UC1  investigation  of  exhaust  pollutant  and  • 
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plume  formation  in  continuous  combustion  is  a combined  analytical 
and  experimental  study  of  turbulent,  backmixed  combustion. 
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Experiments  are  being  conductid  uf.ing  laboratory  backnixed 
combustors  operating  with  premixed  and  nonpremixed  methane/ 
air  and  propane/air  reactants.  Operation  on  liquid  fuels 
is  currently  being  added.  Numerical  models  arc  being  developed 
and  systematically  tested  against  the  experimental  results  to 
evaluate  and  improve  the  models.  The  goal  is  to  identify, 
understand,  and  control  pollutant  and  plume  formation  in 
backnixcd  combustion  typical  oi  gas  turbine  engines. 

The  investigation  is  divided  into  the  following  elements: 
The  implements t ion  of  the  elements  is  detailed  .in  the  Perfornar.ee 
Schedule  presented  in  Figure  2. 

° Model  Development  and  Evaluation 

Evaluate  and  refine  numerical  procedures  and  models 
of  turbulence  by  comparing  numerically  predicted 
profiles  of  velocity,  turbulence  energy,  temperature, 
and  tracer  concentration  of  experimentally  determined 
profiles  (nonreacting  and  reacting  flow). 

Evaluate  and  refine  the  numerical  procedures  and 
the  coupled  models  of  turbulence  and  chemical 
kinetics  by  comparing  numerically  predicted  profiles 
of  velocity,  turbulence  energy,  tempera? ere  * tracer 
concentration,  and  mass  fractions  of  hydrocarbons, 
nitric  oxides,  carbon  monoxide,  oxygen,  and  carbon 
dioxide  to  experimentally  determined  profiles 
(reacting  flow). 

o Mcchetnisas  of  Pollutant /Plume  Formation 

Perform  parametric  studies  (theoretical  and  experi- 
mental) to  identify  the  relative  contribution  of 
the  chemical  reactions,  transpet t processes  (heat 
and  mass  diffusion,  fluid  mechanics,  and  recircula- 
tion behavior),  and  system  cara«aet»rs  (o.g.  geometry 
flew  rates)  to  pollutant  formation. 

o Supplemental  Studies 

Lean  Combustion.  Evaluate  the  emission  character- 
• istics  of  the  opposed  jet  combustor  (OJC/  and  the 

potential  of  utilising  the  OJC  in  stabilizing  lean 
combustion . 
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Swirl.  Design  nrd  Inst  .nil  the  capability  to 
generate  and  conn ol  swirl  in  the  laboratory 
combustor  systems. 

Air  Preheat.  Install  the  capability  of  a.r 
preheat  to  600*C  to  wore  effectively  simulate 
gar.  tuibine  combustion  behavior. 

Non-Preeixed  Combustor.  Expand  the  existing 
premised  combustion  studies  to  include  non-' 
p remixed  systems  norc  representative  of  the 
majority  of  £«s  turbine  systems. 

Par  t I cul  at  o Forga  t ion  /l.unlnosl  t y . Investigate 
particulate  formation  in  reacting  flows  with 
recirculation  with  special  emphasis  on 
alternative  fuels.  Investigate  luminosity  in 
reacting  flows  with  recirculation  with  special 
attention  to  alternative  fuels  and  inpact  on 
combustion  liner  degradation. 

NO  Probe.  Explore  chemical  transf ornat ion  of 
nitrogen  oxides  in  sample  probes  and  sample 
lines  while  sampling  combustion  products. 

Such  information  is  important  in  the  use  of 
measured  concentration  data  in  pollutant  and 
plume  formation  studies. 

3.4  The  DC  I Combustor  3 ys teas 

Two  combustor  systems — premixed  and  nonpre.uixec — arc 
employed  in  the  UCI  wtudy.  Schematics  of  both  systems  arc 
presented  in  Figure  3.  The  premixed  system,  shown  in  Figure 
3a,  is  an  opposed-reacting- jet  combustor  (OdC) . The  non- 
prer.ix*»d  system,  shown  in  Figure  3b,  is  a classical  con- 
figuration with  fuel  injected  at  the  axis  in  a stream 
concurrent  with  the  surrounding  air. 

The  opposed-jet  combustor  (OJC)  exhibits  essential 
features  found  in  practical  continuous  combustion  systems 
(backmixing,  high  intensity  combustion)  but  avoids  the  attendant 
coeplev.i ties  (nonpremixed,  two  phase  combustion)  appropriate 
for  complementary  studies.  The  recirculation  none  is  isolated 
and  free  from  geometrically  governed  boundary  conditions. 
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The  isolation  .allows  the  tu rbu lent /hydrodynamic  effe  t s to 
be  explored  apart  froe  the  complicating  questions  associated 
with  boundary  condition  specification.  The  features  of  the 
OJC  especially  attractive  and  novel  to  the  study  are  the 
extended  range  and  control  of  turbulence  intensity,  nixing 
lengths,  and  nixture  equivalence  ratios  versus  those  offered 
by  the  sudden  expansion  on  physical,  bluff  body  configurations. 
For  cxiar-lc,  the  rate  and  intensity  of  backnixing  in  the  re- 
circulation zone  nay  be  altered  experimental ly  and  numerically 
over  a broad  range  by  varying  the  ratio  of  approach  velocity 
to  let  velocity. 

The  CiC  combustor  i*  a 51 mm  I.S.  x 457nn  cylindrical  Vycor 
chamber  containing  ac  aerodynamic  (opposed- jet)  f X3me-holder . 

The  incoming  flow  cf  prefixed  fuel  (methane  nr  prepane}  is 
opposed  by  a high  velocity  jet  (£,/»<<  5)  i*s»iing  from  a 

j a 

1.3ma  1.8.  (G . 4cm  O.D.)  up tcr-cof/led  tube.  The  jet  is  coin" 
cident  with  the  cocbuotor  axis  and  located  7 upstrena  iroa 
the  combustor  exit.  The  backmixed  cone  enquired  for  flame 
stabilization  is  generated  along  the  jet  boundary,  The  cota- 
bustor  is  currently  bvting  enlarged  to  ?Sma  to  ieptove  the 
resolution  of  the  flovfinld  spatially. 

The  nonpreaiixed  system  Is  currently  being  designed-' 

The  combustor  size  will  match  the  expanded  pramlxcd  version 
of  78na.  The  nonpremixed  configuration  viJl  allow  exploration 
of  the  effects  of  mixing  the  two  dissimilar  streams  on 
pollutant  formation  and  numerical  prediction,  allow  utilization 
of  liquid  fuel,  end  provide  a bridge  between  the  premixed  system 
and  the  more  conventional  turbine  combustor  configurations. 
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2.  BACKGROUND 
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2.1  Numerical 

The  devclopnent  of  a predictive  nethod  for  continuous 
conbustion  must  be  founded  upon  a suitable  conputat ional  pro- 
cedure. Because  the  flows  are  characterized  by  strong  back- 
mixing,  the  describing  partial  differential  equations  must  be 
elliptic. 

A computational  technique  nay  be  developed  to  solve  the 
elliptic  equations  wherein  the  goal  is  to  solve  directly  for 
the  instantaneous  velocity.  This  approach  is  costly, 
relatively  conplcx,  and  requires  substantial  conputer  storage 
capacity.  An  alternative  approach  for  engineering  calculations 
is  to  adopt  tine-averated  properties.  Here,  Instantaneous 
properties  (pressure,  tenperaturc,  velocity.,..)  arc  replaced 
py  the  sun  of  the  mean  (tine-averaged)  value  of  the  property 
and  a fluctuating  component.  The  governing  equations  arc 
then  structured  in  terns  of  the  tine-averaged  properties  plus 
additional  terns  involving  the  fluctuating  components.  The 
problca,  then,  is  transforned  from  solving  for  the  tine- 
dependent  instantaneous  velocity  to  providing  a codel  that 
account*  for  the  nev  fluctuating  conponcnt  terns  in  the  tiae- 
averaged  equations. 

One  Approach  used  to  account  for  the  fluctuating  component 
terns  i<*  to  oodify  the  transport  properties  viscosity, 

ji)  . Thu  laodificd  values  of  the  transport  properdins  ere  called 
Vf*  properties  (e . . effective  viscosity,  on<* 
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formulated  to  account  for  the  turbulence  transport  associated 
with  the  fluctuating  components.  Examples  of  two  effective 
viscosity  models  are  the  following  algebraic  and  two-equation 
formulations : 


i 
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Algebraic3 


Pcf,  - K D*/J  L_1/s  pJ/s  (nV2)1/1  (1) 

K:  empirical  coefficient 

D:  diameter  of  combustor 

L:  length  of  combustor 

p:  local/ density 

mV2:  inlet  kinetic  energy 


Two-Equation  K 

“off  * cp  P k’'c  (2) 

c^:  empirical  coefficient 

p : local  density 

k : local  turbulent  kinetic  energy 

C : local  turbulent  energy  dissipation 

rate 

For  a given  geometry  and  inlet  flow  rate,  the  algebraic 
formulation  is  functionally  dependent  only  upon  local  density. 
The  two-equation  model  introduces  additional  unknowns  that 
complicate  the  numerics  and  adds  computational  cost  penalties. 
However,  the  tvo-cquation  effective  viscosity  model  is 


m 


t 


I 
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intuitively  aorc  accurate  because  it  accounts  for  spatial  * 
variations  in  turbulent  kinetic  energy  and  length  scale. 

The  turbulent  transport  of  scalar  quantities  is  related 
to  the  eddy  viscosity  through  the  appropriate  Prantitl  or  Scheldt 
numbers  as  follows: 


vef  f 
Ceff,* 


reff,<>! 
$ : 


effective  exchange  coefficient  for 
dependent  variable  $ 

Prandtl  or  Schmidt  number 
dependent  variable 


(3) 


Table  I lists  the  following  major  computer  codes  that  have 
been  developed  using  the  tine-averaged  approach  to  solve  the 
elliptic  flow  equations: 

PISTEP  II  is.  a </-u>  progran  with  a_  two-equation 
effective  viscosity  model.  The  UCI  version  includes 
cosbusticn  chcnistry,  the  OJC  geometry,  and  the  algebraic 
description  of  effective  viscosity  for  convenience 
in  initiating  a first  solution. 

TEACH  is  a three-dimensional  program  written 
with  p-v  as  the  major,  dependent  variables.  The 
code  has  the  attractive  features  of  a third  dimen- 
slop  option  and  an  efficient  iterative  routine. 
e The  UCI  version  includes  both  a two**paranetcr  and 

algebraic  description  of  the  eddy  viscosity  models 
as  well  as  the  combustion  chcnistry,  and  OJC 
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.geometry  — 

A third  code,  CRISTY,  developed  at  United 

Technology  Research  Center  (UTRC)  through  a 
Joint  sponsorship  by  the  Environmental  Protection 
Agency  (EPA)  and  the  Air  Force  Aero  Propulsion  Lab- 
oratory (AFAPL) , is  being  run  in  cooperation  with 
the  ERA. 

2.2  Experimental 

The  experimental  dimension  of  the  UCI  program  combines 
the  following  components  of  the  continuous  flow  combustion  test 
facility  in  the  UCI  Combustion  Laboratory: 

• combustor  test  facility 

• gas  analysis  system 

• laser  veloclmeter 

The  integrated  system  represents  the  unique  experimental 
capability  of  providing  detailed  measurements  of  composition, 
temperature,  velocity  and  turbulence  data  to  (1)  develop  and 
evaluate  the  predictive  methods,’ (2)  explore  the  mechanisms  of 
exhaust  pollutant  and  plume  formation,  and  (3'  conduct  supplemental 
studies. 

Combustor  Test  Facility.  The  test  facility  includes 
built-in  controls  and  measuring  devices  for  monitoring  operating 
parameters  as  veil  as  diagnostics  for  the  flame  and  emission 
characteristics.  Provisions  are  included  for  varying 
operational  parameters  continuously  while  the  combustor  is 
firing.  The  structural  features  of  the  combustor  test  section 
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readily  accommodate  mechanical  probing  for  sample  input  to 
the  gas  analysis  system  and  optical  probing  by  the  laser 
velocimeter.  Combustor  sir  inlet  temperatures  up  to  600*C  . ■ 

(1112*F)  can  be  accommodated.  A flow  diagram  of  the  test 
facility  is  sliovn  in  Figure  4.  The  combustor  test  section 
is  nominally  51mm  in  diameter  and  is  currently  being  expanded 
to  78mm.  The  air  is  supplied  on  a continuous  basis  (O.lKg/sec) 
vith  approach  velocities  variable  to  18,3  netcrs/sec  (60  ft/sec). 


The  combustor  facility  is  interfaced  with  .an  automated 
data  acquisition  system  coupling  an  analog  data  scanner  with 
a minicomputer.  System  capabilities  include  real-time 
monitoring  and  control  of  the  experiments  and  efficient  data 
reduction.  A block  diagram  of  the  data  acquisition  system' 
is  shown  in  Figure  S. 

Gas  Analysis  System.  A packaged  turbine  engine  analysis 

system  (Scott  Research  Labs.,  Model  113)  is  available  for  the 

exhaust  gas  analysis.  The  unit  consists  of  five  commercially 

available  instruments:  Two  nondispersivc  infrared  analyzers 

(NDIR)  for  measurement  of  CO  and  COj,  a chemi) ucinesccnt 

analyser  (ChA)  for  measurement  of  HO/NO^,  a flame  ionization 

detector  (FID)  for  total  hydrocarbons,  and  a paramagnetic 

oxygen  analyzer.  Also  included  is  the  associated  peripheral 

equipment  necessary  for  sample  conditioning  and  instrument 

calibration.  All  instruments  supply  an  analog  output  to  a 

5-channcl  strip  chart  recorder  and/or  the  data  acquisition 
* 

systems.  A detailed  breakdown  of  the  analytical  instruments 
is  shown  in  Table  II. 


Laser  Vclocloeter.  A DISA  type  SSL  laser  Doppler  veloci- 
neter  (LDV)  system  Is  available  for  velocity  measurement  In  the 
combustor  flovfield.  The  principal  components  of  the  LDV  system 
include  a 5 nW  Hc-Ne  laser  light  source  manufactured  by  Spectra- 
Physics  (Model  120)  an  integrated  optical  unit  and  a DISA 
(Model  55L20)  frequency  tracker  for  Doppler  signal  processing. 
The  signal  processing  is  currently  being  expanded  to  include 
frequency  counting  through  use  of  the  PDP-11  computational 
and  data  sampling  capability. 

The  LDV  optical  bench  is  mounted  on  a motorized  traverse 
to  facilitate  rapid  scanning  of  the  flovfield  while  maintaining 
an  accurate  spatial  resolution.  A schematic  diagram  of  the  LDV 
optical  arrangement  and  associated  electronics  is  shown  in 
Figure  6. 

In  our  experience  to  date,  it  has  been  found  necessary 
to  artifically  seed  the  flow  system  with  supplementary  light 
scattering  particles  to  provide  an  adequate  signal  input  to 
the  data  processing  instrumentation.  The  nonreacting  flow 
studies  arc  conducted  using  salt  (NaCl)  particles  introduced 
into  the  flow  as  an  aerosol.  Combustion  experiments  necessitate 
introducing  high  molting  point  solids,  such  as  alumina  A^O-j, 
by  a fluidized  bed  technique. 

2.3  Coupled  Huocrlcol/Expcrimcntnl  Procedure 

The  program'  plan  of  the  investigation  is  designed  to 
accomplish  the  objectives  of  the  study  via  a judicious  coupling 
of  numerical  methods  to  experiment.  The  approach  is  to  address 
in  sequence  isothermal  (nonreacting),  heated  (nonreacting  with 
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preheat),  and  hot  (reacting  flow),  The  approach  ia  outlined 
in  Table  II  and  illustrated  in  Figure  2.  The  eddy  viscosity 
and  mass  transport  submodels,  and  initial  conditions  and 
boundary  conditions  arc  first  tested  against  experiment  for 
isothermal  flow.  The  tests  used  include  comparisons  of  vel- 
ocity profiles,  turbulence  intensity  profiles,  and  profiles 
of  tracer  concentrations.  The  eddy  viscosity,  mass  transport, 
and  energy  transport  subaodcls,  and  initial  and  boundary 
conditions,  arc  secondly  tested  against  experinent  for 
heated  (nonreacting)  flow.  The  main  stream  and-  jet  stream  are 
independently  heated  to  temperatures  in  excess  of  400*C  (754*F). 
The  tests  used  include  comparisons  of  velocity  profiles,  tur- 
bulence Intensity  profiles,  profiles  of  tracer  concentration, 
and  temperature  profiles.  Third,  the  eddy  viscosity,  mass 
transport,  energy  transport,  and  chemistry  submodels,  and  initial 
and  boundary  conditions  are  tested  against  experiment  for  hot 
(reacting)  flow.  The  tests  include  those  above  plus  comparison 
of  species  concentration  profiles  of  the  major  products  of  com- 
bustion as  well  as  carbon  monoxide  and  oxides  of  nit'ogen. 

The  rate  and  intensity  of  backmixing  in  the  r>  -culating 
flow  region  are  altered  numerically  and  experiments ■ Ij""by 
variations  of  the  ratio  of  approach  velocity  and  jet  velocity. 
Turbulence  intensity  and  scale  (eddy  size)  are  altered  by  the  use 
of  screen  grids  upstream  of  the  stabilization  zone.  The  com- 
parative tests  encompass  a broad  range  of  turbulent  transport 
rates  between  the  recirculation  zone  and  the  primary  gas  stream. 


Additional  parametric  variations  include  the  mixture  equivalence 
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2,4  HO  Probe 


Measurement  of  combustion  products  in  laboratory  burners 
and  flames  is  required  to  identify  the  elements  of  chenical 
kinetic  and  fluid  mechanics  responsible  for  the  formation 
of  pollutant  species. 

Measurement  of  the  gaseous  product  composition  emitted 
from  a combustion  source  generally  proceeds  by  extracting  an 
exhaust  sample  for  subsequent  quantitative  analysis.  In  order 
to  obtain  reliable  data,  precaution  must  be  taken  to  ensure 
that  the  analytical  instrumentation  receives  a sample  that 
is  chemically  identical  to  the  composition  existing  at  the 
sampling  point.  Potential  sample  transformations  may  be 
minimized  by  careful  selection  of  materials  contacting  sample 
gases  in  the  physical  probe  used  for  sample  removal  and  in 
the  line  used  to  transport  the  sample  to  the  instrumentation. 

Of  all  the  significant  products  of  combustion,  nitrogen 
oxides  (NO^)  are  especially  suceptlble  to  chemical  transfor- 
mation. In  addition  to  the  total  nitrogen  oxides  concentration, 
the  concentration  of  the  primary  species  in  the  nitrogen  oxides 
family,  nitric  oxide  (HO)  and  nitrogen  dioxide  (NO^),  must 
often  be  determined.  For  example,  spatial  distributions  of 
nitric  oxide  (HO)  and  nitrogen  dioxide  (N^)  are  re<lu*re<*  to 
identify  the  chemical  kinetic  mechanisms  responsible  for  the 
formation  of  nitrogen  oxides  in  flames  or  gas  turbine  combustors. 
The  oxidation  of  HO  to  NO^  or  the  reduction  of  NO^  to  HO  in 


\ 
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sample  probea  anil  sample  lines  prevents  an  accurate  determination 
of  eitlier  110  or  HC>2.  As  a result,  oxidation  or  reduction 
reactions  involving  NO  or  N02  may  often  be  unacceptable  even 
though  total  NO^  is  conserved. 

An  adequate  accounting  of  NO^  transformations  requires 
that  experiments  be  conducted  to  identify  (1)  the  chemical 
transformations,  if  any,  that  occur,  (2)  the  extent  to  which 
they  occur  during  typical  sampling  times,  and  (3)  the  conditions 
(e.g.  sample  temperature,  sample  composition,  sample  line 
material)  for  which  they  occur. 

The  present  experimental  study  is  designed-  to  assess  NO^ 
transformations  that  may  be  encountered  when  sampling  within 
flames  and  at  the  exhaust  plume  of  laboratory  combustors.  The 
experimental  system  is  currently  designed  to  simulate  the  actual 
conditions  experienced  in  sample  lines  used  for  sampling  gaseous 
combustion  products.  A diagram  of  the  experimental  system  is 
shown  in  Figure  7.  Test  parameters  include  composition  of 
the  carrier  gas,  dopant  gases,  and  probe  material. 

A carrier  gas  simulating  the  primary  combustion  products 

is  selected  from  one  of  three  prepared  sources  of  0,  1,  or 

SZ  02>  12Z  C02,  and  balance  Nj.  The  carrier  gas  flow, 

4 liters/minute,  is  doped,  with  CO,'  NO,  and  NO,,  metered  from 

high  concentration  source  cylinders  by  means  of  porous  sintered 

metal  flow  restrictors.  After  doping,  the  carrier  gas  enters 

a silica  preheat  tube  that  raises  the  gas  temperature  to  the 

desired  probe  best  temperature.  From  this  point,.,  the  doped 

carrier  gas  enters  a sample  probe  test  section.  An  oven  is 
♦ 

used  to  maintain  the  temperature  of  the  doped  carrier  goo  at 
the  desired  probe  test  temperature. 


li 
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The  temperature  of  the  gas  vithin  the  sample  probe  Is 
increncntally  varied  fron  25®C  to  400*C,  Tcnpcraturcs  of  the 
gas  stream  (T^  and  T^)  arc  measured  with  insulated  platinum 
resistance  thermometers  centered  in  the  probe  bore  at  the  it\let 
and  outlet  of  the  sample  probe.  The  oven  temperature  is  also 
recorded  by  a thermocouple  located  adjacent  to  the  outer  dia- 
meter of  the  sample  probe.  The  sample  probe  heater  can 
accommodate  a variety  of  probe  materials  of  lengths  up  to  2 


Cas  composition  is  measured  at  sample  Point  3 to  assess 
the  extent  of  NO  and  HO^  transformation  within  the  2 meter 
sample  probe  segment.  Sample  lines  leading  from  points  2 
and  3 (Figure  7)  are  short,  equal-distant,  and  made  of  1/4 
inch  o.d.  TFF.  Teflon.  Screening  tests  using  varying  lengths 
of  TFE  Teflon  have  been  conducted  to  assure  that  NO^  absorption 
is  not  a significant  factor  in  the  present  experiment. 

Analysis  of  NO  and  NO^  is  conducted  with  a Beckman  Model 
951H  chemiluminescent  oxides  of  nitrogen  analyzer,  NO^  is 
determined  by  the  difference  between  the  measured  NO  and  NO 
concentrations.  The  (carbon)  converter  efficiency  is  moni- 
tored by  periodic  tests  using  the  methods  outlined  in  the 
federal  Register.5  Measurements  of  carbon  monoxide  concentration 
is  conducted  v*th  a Beckman  Model  315BL  nondlspcrsive  infrared 


analyzer. 
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3.  DESCRIPTION  OF  PROGRESS : REPORTING  PERIOD 


The  performance  schedule  is  presented  in  Figure  2.  The 
first  ten  months  of  the  research  program  were  structured  to 
accommodate  (1)  modification  and  testing  of  the  experimental 
system,  (2)  modification  of  the  numerical  methods  to  the  geo- 
metry of  the  opposed  reacting  combustor,  (3)  expansion  of  the 
numerical  methods  to  include  multicomponent  reacting  flow  and 
(4)  initiation  of  the  cold  flow  numerical  and  experimental 


The  current  reporting  period  covers  the  subsequent 
twelve  (12)  month  period  and  Includes  (1)  the  test  of  the 
numerical  methods  against  experiment  for  the  case  of  isothermal 
flow,  (2)  expansion  of  the  numerical  methods  to  include  propane 
oxidation  kinetics,  (3)  expansion  of  the  experimental  facility 
to  Include  air  preheat  capability,  (fi)  initiation  of  the  gen- 
eration of  a broad,  experimental  data  base  for  the  cases  of 
isothermal  and  hot  flow,  and  (5)  initiation  t»f  an  in-depth 
exploration  of  chemical  transformation  that  may  occur  in  probes 
and  lines  used  to  sample-  and  transport  nitrogen  oxides. 

The  program  is  on  schedule. 

3.1  Numerical 

A summary  of  the  improvements  and  modifications  completed 
during  the  current  reporting  period  is  outlined  in  Table  IV. 

^Current  problems  outstanding  and  plans  for  their  resolution  arc 
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also  identified. 
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PISTEP  II  and  TEACH  were  tested  against  experiment  for 
the  case  of  isothermal  flow.  The  results  of  the  TEACH  eval-’ 
uatlon  are  discussed  in  Section  3-3  (Combined  Numerical/ 
Experimental)  and  will  be  presented  at  the  Sixtecntl  Symposium 
(International)  on  Combustion.  The  PISTEP  II  evaluation  is 
currently  being  completed. 

In  cooperation  with  the  EPA,  the  CRISTY  Code  was  modified 
to  the  UCI  OJC  geometry  and  is  currently  being  -run  for  the 
algebraic  eddy  viscosity  model. 

Solutions  have  been  obtained  for  hot  (reacting  flow) 
mixtures  of  methane/air  with  nitric  oxide  formation  using 
TEACH.  Solutions  are  now  being  pursued  with  PISTEP  II.  Both 
codes  arc  undergoing  modifications  to  consider  propane  oxidation 
kinetics. 

In  summary,  PISTEP  IX,  TEACU,  and  CP.ISTY  are  specialized 
to  the  OJC  geometry , and  are  now  being  tested  for  the  case  of 
isothermal  conditions.  The  evaluation  of  the  methods  for  the 
ease  of  hot  flow  has  comsensod.  The  codes  arc  being  expanded 
to  include  propane  kinetics  as  veil  as  methane  oxidation 
chemistry.  Tests  of  the  codas  for  the  hot  flow  condition  will 
continue  during  the  continuation  year. 

3. 2 Experiment 

Combustor  Test  Facility.  A summary  of  all  the  improvements 
-and  modifications  completed  during  the  current  reporting  period 
are  outlined  in  Table  V.  Current  problems  outstanding  and  plans 
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for  their  resolution  are  also  identified. 


The  major  modification  to  the  combustor  test  facility'  ' 
was  the  inclusion  of  air  preheat  capability  to  600*C  (1112*1'). 


Electric  heating  was  selected  (1)  to  avoid  contamination  of 


the  coabusior  inlet  air  (contamination  of  the  inlet  air  could 


significantly  alter  the  combustion  process  chemistry),  (2) 


to  provide  a precise  control  of  tcmpcrscui e,  and  (3)  to  allow 


a rapid  response  to  changing  flow  conditions. 


The  combustor  air  preheat  system  is  designed  to  have  the 


flexibility  to  heat  the  combustoi  inlet  air  te  any  temperature 


virhin  the  design  temperature  range  and  to  be  operational  for 


a range  of  airflows.  The  maximum  design  temperature  for  the  air 


preheat  system  is  600°C.  Selection  of  this  temperature  is 


based  on  actual  gas  turbine  operating  conditions.  The  air  preheat 


system  will  provide  heated  inlet  air  fur  combustor  air  velocities 


c£  7.62  to  15.24  m/scc  (25  to  50  feet/sec). 


The  Initial  scries  of  combustor  preheat  studies  will  involve 


the  preoixed  combustion  of  propane  With  heated  air.  The  ignition 


temperature  for  a f uel-cxidanr.  mixture  is  the  temperature  at 


which  self  sustained  combustion  will  occur.  Reported  ignition 


temperatures  for  propane-air  mixtures  range  from  450"C  to  600*C  . 


To  avoid  the  possible  hazardous  condition  of  mixing  propane  with 


air  heated  to  the  ignition  temperature,  a maximum 
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operating  temperature  of  A00*C  will  be  observed.  This  maximum 
temperature  will  allow  combustor  studies  to  be  performed  at 
temperatures  typical  oi  low  pressure-ratio  gas  turbine  engines. 
Future  combustion  studies  using  a nonpreraixed  combustor  con- 
figuration will  allow  the  maximum  combustor  air  preheat  system 
temperature  to  be  extended  to  the  desired  limit  of  600*C. 

Laser  Vclociacter,  The  LDV  development  has  emphasized 
isothermal  flow  conditions,  A commercially  available  nebulizer 
is  used  to  introduce  supplementary  laser  light  scattering 
particles  into  the  flow.  The  nebulizer  atomizes  high  concen- 
trations of  uniform  size,  salt-water  solution  droplets.  The 
aerosol  is  dried  to  small,  solid  salt  particles  of  uniform  size 
before  entering  the  combustor.  The  size  of  the  particles  is 
controlled  by  the  original  salt  concentration  in  the  solution. 

Current  problems  outstanding  concern  signal  acquisition 
in  the  recirculation  zone  and  the  stress  layer  located  along 
the  peripherty  of  the  opposed  jet.  Several  approaches  are 
currently  being  investigated  to  improve  data  acquisition  and 
interpretation.  Bnphasis  is  being  placed  on  development  of  a 
period  timed  (counter)  signal  acquisition  system  to  complement 
the  current  frequency  tracking  technique. 

3.3  Coupled  Expcrlnental/Kuacrlcnl 

One  ot  the  important  elements  of  the  current  research 
program  is  the  coupling  of  the  numerical  and  experimental 
components.  The  development  of  a realistic  numerical  model  of 
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continuous  combustion  requires  careful  validation  by  vs y of 
cxpcrincntdl  observation. 

Emphasis  during  the  current  reporting  period  baa  bren 
directed  to  (1)  testing  the  numerical  m*\.hods  for  the  case  of 
isothermal  flow  (to  evaluate  the  eddy  viscosity  aod<?l)  and, 

(2)  testing  the  nuacrical  methods  for  the  case  of  isothermal 
flow  with  a tracer.  For  the  case  of  isothermal  flav,  the 
eddy  viscosity  submodels  were  first  evaluated  for  the  prediction 
of  system  hydrodynamics.  A tracer  gas  was  later  introduced  to 
evaluate  the  eddy  viscosity  submodels  for  nonuniform,  multi- 
component flow  where  mass  as  veil  as  momentum  exchange  is 
Important.  Experimental  measurements  and  numerical  predictions 
of  velocity  and  tracex.  isopleths  for  a range  of  approach 
velocities  and  backmix  characteristics  provided  the  critical 


The  results  reported  here  emphasize  the  ease  of  isothermal 
flow.  A satisfactory  assessment  and  understanding  of  the  iso- 
thermal case  is  a prerequisite  to  engaging  the  hot  flow  ease. 
Combustion  Imposes  the  additional  requirement  for  submodels  of 
chemistry  and  turbulence/chcm.lstry  interaction  that  may  mask 
fundamental  deficiencies  in  the  eddy  viscosity  submodel. 

However,  hot  flow  results  are  presented  as  an  indication  of  the 
suitability  of  the  two  eddy  viscosity  submodels  in  the  pre- 
diction of  reacting  flows. 

The  isothermal  and  hot  flow  results  collectively  demonstrate 
the  merits  of  two  eddy  viscosity  submodels  for  the  prediction 
of  transport  behavior  in  continuous  combustion  systems  of 
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practical  engineering  interest  (turbulent  reacting  flows  *j 

with  recirculation).  In  addition,  the  results  provide  useful  j 

information  for  the  development  of  applied  computational  topis 
for  analyzing  the  performance  of  continuous  combustion  systems. 

The  present  study  explores  specification  of  the  eddy 
viscosity,  turbulent  mass  transport,  in  con- 

i ^ 

fined  flows  exhibiting  strongbackmixing  or  recirculation. 

Comparative  calculations  were  performed  using  two  distinct 
turbulence  submodels  to  define  the  eddy  viscosity — an  algebraic 
model  (Equation  1)  2nd  a more  elegant  two-equation  submodel 
(Equation  2) . 


RESULTS  (ISOTHERMAL  FLOW) 

To  effectively  test  the  submodels  of  eddy  viscosity 
over  a range  of  turbulent  conditions,  experimental  data  were 
obtained  for  approach  velocities  of  15,24  and  7.62  m/sec  and 
jet  velocities  of  30,5,  61  and  130n/sec,  Velocity  profiles 
were  mapped  throughout  the  flovfield  using  a 1mm  0,D.  pitot 
tube.  Tracer  profiles  were  generated  by  introducing  altern- 
atively pure  carbon  monoxide  (CO) 'in  the  jet  and  pure  carbon 
dioxide  (CO^)  in  the.  mainstream.  The  mixing  between  the  jet 
and  mainstream  was  assessed  by  mapping  CO  and  CO2  profiles 
throughout  the  flovfield  using  a 3am  stainless  steel  tube 
and  Beckman  3156  ND1R  instruments  for  the  gas  sampling  and 

analysis.  The  Reynolds  number  of  the  approach  stream  ranged 

* 4 4 

from  2.5  x 10  (at  7.62n/sec)  to  5.0  x 10  (at  15.24n/sec). 
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tiomcntion  Transport 

The  eddy  viscosity  submodels  were  first  tested  for  a homogeneous, 
isothcrnal  flowficld  to  explore  the  effectiveness  of  the  submodels  in 
predicting  momentum  transport.  Experimental  and  predicted  profiles  of  the 
axial  component  of  mean  velocity  are  presented  in  Figure  8 for  an  average  inlet 
velocity  of  15.24n/scc.  The  selected  profiles  emphasize  the  flow  regime 
dominated  by  recirculation.  The  calculated  velocity  profiles  for  both 
turbulence  submodels  assume  an  experimentally  determined  velocity  profile 
as  an  inlet  condition. 

The  trend  of  the  velocity  profiles  agrees  qualitatively  with  the 
experimental  results.  The  two-equation  submodel  retains  the  Inlet  turbulent 
velocity  profile  in  the  approach  section  while  the  velocity  profiles 
predicted  by  the  algebraic  submodel  flatten  out  due  to  the  spatially 
uniform  viscosity.  The  presence  of  the  opposing-jet  is  communicated  further 
upstream  for  the  tvo-equaticn  submodel  than  for  the  algebraic  eddy  viscosity, 
and  the  two-equation  submodel  demonstrates  better  correlation  with  the 
experimentally  observed  stagnation  point. 

The  mean  velocity  profiles  in  the  rccirculaticn  zone  are  also  better 
represented  by  the  two-cquation  turbulence  submodel,  although  both  sub- 
models over-predict  the  jet  expansion.  In  addition,  the  two-cquation 
velocity  profiles  decrease  too  rapidly  near  the  chamber  wall. 

Downstream  from  the  backmixed  zone  both  submodels  fail  to  adjust  for  the 
redistribution  of  the  bulk  flow  in  the  turbulent  wake  of  the  recirculation 
zone  and  under-predict  the  velocity  adjacent  to  the  Jet  wall.  Both  sub- 
models display  good  agreement  in  the  bulk  flow  regime  in  the  absence  of 


recirculation  or  wall  effects. 


24 


F " ' 


! 

I 

4 


j 

I 

t 

i 


i 


j 


i 

i 

\ 


t 


i 


The  recirculation  zone  strength  was  examined  by  reducing  the  approach 
velocity  by  a factor  of  two  while  retaining  the  sane  jet  velocity.  The 
cxpcrlncntal  and  predicted  velocity  profiles  arc  presented  in  Figure  9.  . 

The  calculated  results  demonstrate  a significant  departure  from  experiment 
in  the  region  of  strong  backmixing.  The  reverse  flow  region  again  penetrates 
further  upstreaa  in  the  two-equation  submodel  prediction. 

The  acceleration  of  the  bulk  flow  around  the  recirculation  zone  is 
under-predicted  in  contrast  to  the  over-predicted  jet  expansion  for  both 
submodels.  This  result  is  especially  evident  near  the  plane  of  the  jet 
exit  where  the  redistribution  of  the  bulk  flow  with  the  jet  discharge  is 
restricted.  This  effect  is  conveyed  downstream  from  the  recirculation 
zone  where  the  tvo-cquatlon  velocity  profiles  correlate  well  with 
experiment  near  the  chamber  wall  but  deviate  near  the  jet  wall. 

Hass  Transport 

The  tve  eddy  viscosity  submodels  were  secondly  tested  in  conjunction 
with  Che  turbulent  mass  exchange  submodel  (Eq.  3)  for  the  conditions  of  an 
isothermal,  nonhomogeneous  flowfield.  Spatial  distributions  of  the  predicted 
tracer  concentrations  are  compared  to  experimental  measurements  in  Figures  10 
and  11  for  the  baseline  conditions  outlined  previously. 

Figure  10.  presents  the  algebraic  and  "two-equation  submodel  predictions, 
and  experimental  tracer  concentration  distributions  for  a 15.24m/scc 
average  inlet  velocity.  (The  jet  stream  is  100  percent  carbon  monoxide.) 
Immediately  evident  is  the  contrasting  axial  and  radial  spread  of  the 
predicted  tracer  concentration  profiles.  Differences  in  axial  transport  of 
CO  arc  consistent  with  the  observations  forwarded  in  the  previous  section 
on  noncntua  transport.  Hore  notable  arc  the  diverse  radial  transport  pro- 
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dictions.  The  spreading  of  the  jet  discharge  into  the  mainstream  for  the 
algebraic  eddy  viscosity  conforms  to  tlie  experimental  results. 

In  contrast,  the  two-equation  submodel  Units  the  radial  nixing  of  the 
jet  atreara  (tracer)  with  the  bulk  flow.  The  segregation  of  the  jet  and 
the  mainstream  is  verified  by  the  high  CO  concentrations  along  the  jet 
vail,  in  contrast  to  the  absence  of  CO  near  the  chamber  wall.  This  result 
was  further  substantiated  by  introducing  a carbon  dioxide  (C02)  trorcr  in 
the  approach  stream.  The  numerical  analysis  of  CO ^ concentrations 
for  the  two-equation  submodel  predict  negligible  mass  transport  of 
(X>2  from  the  mainstream  into  the  recirculation  2one.  Downstream,  the 
predicted  mixing  of  C02  from  the  bulk  stream  into  the  flow  associated 

with  the  jet  stream  is  minor.  Experimentally  C02  was  observed  to  be 
well  mixed  into  the  jet. 

The  results  for  the  strong  recirculation  zone  .(Figure  11)  demonstrate 
similar  trends.  Although  the  velocity  data  presented  earlier  in  Figure  9 
gave  poorer  agreement  for  conditions  of  strong  recirculation,  the  mixing 
is  improved  because  of  the  enhanced  backnixing. 


RESULTS  (HOT  FLOW) 


A two-step  global  reaction  mechanism  for  methane  oxidation  was 
adopted  for  the  present  study  to  illustrate  the  applicability  of  current 
numerical  methods  to  the  prediction  of  continuous  combustion  flows, 
and  to  briefly  explore  the  suitability  of  the  eddy  viscosity  submodels 
for  the  case  of  hot  flow: 


a,4  + °2 


CO  + HjO 


CO  + 0, 


A detailed  description  of  the  numerical  formulation  of  the  reaction  rate 

expressions  and  boundary  conditions  may  be  found  elsewhere*.  The  dependent. 

variables  for  the  hydrocarbon  system  include  the  mass  fractions  Y_,  and 

Y^  . Distributions  of  the  other  major  species-r-water  (H^O),  oxygen  (op, 

carbon  monoxide  (CO) — are  related  to  Y_„  and  Y by  elemental  mass 

4 t"'"2 

conservation.  An  initial  solution  for  the  reacting  flowfield  was  obtained 
using  published  reaction  rates’  which  were  later  refined  to  more  closely 
resemble  the  experimental  data.  This  approach  was  justified  on  the  basis 
that  the  reaction  rate  data  were  obtained  from  the  study  of  a well-stirred  reactor 
and  may  not  apply  to  the  current  study,  which  deviates  substantially  from 
well-stirred  conditions.  Global  reaction  rates  were  judged  acceptable  for 
the  present  study  because  of  the  uncertainties  encompassing  the  specification 
of  the  systen  aerodynamics. 

Since  the  ala  of  the  study  is  to  predict  the  performance  of  backmixcd 
combustion  systems,  the  hot  flow  calculations  were  extended  to  include  the 


prediction  of  the  pollutant  species,  nitric  oxide  (NO).  Nitric  oxide 
kinetics  were  based  on  the  familiar  Zeldovich*  mechanism:  ■ 


0 + N2  — * HO  + N 


N + 02  — • N0..+  0 


Noting  that  reaction  (6)  is  the  rate  limiting  step  and  adopting  the 
simplifying  assumption  of  0/0 2 equilibrium  results  in  the  reaction  rate 
expression: 

j tHQl  . ?k  K fN  1 10  1** 

. dt  2 of  0 *N2Jl02J 


The  temperature  and  02  distribution  obtained  from  the  solution  of  the 


saafldw; 


hydrocarbon  system  was  used  as  a basis  for  the  HO  kinetic  calculations. 

The  experimental  tests  conducted  to  complement  the  numerical  predictions 
of  hot  flow  properties  utilized  stoichiometric  proportions  of  premixed 
methane  and  air  in  both  the  main  and  Jet  streams.  The  reactants  were 
Initially  at  ambient  temperature  and  the  combustion  was  completed  at  atmos- 
pheric pressure.  The  test  matrix  included  the  same  approach  velocities 
selected  ior  the  isothermal  cases  — 7-62  and  15.2Wsec  and  a Jet 
velocity  of  130a/sec.  Detailed  maps  of  flovfield  properties  were  obtained 
by  conventional  probing  techniques.  Temperature  measurements  were  made  using 
an  uncoatcd  Pt/Pt-132  Rh  thermocouple  and  a Scott  Model  12'5  chemiluminescent 
analyzer  was  used  for  oxides  of  nitrogen  measurements.  Gas  samples  were 
extracted  via  a novcablc  3.2mm  O.D.  watercoolcd.  316  stainless  steel  probe 
having  a tubular  inlet,  and  conveyed  through  a heated  teflon  sample  line 
to  a packaged  exhaust  gas  analysis  system.  Selected  results  of  hot  flow 
properties  are  presented  in  Figures  12  and  13  for  the  conditions  indicated. 
Temperature  and  SO  profiles  are  presented  as  indicators  of  useful  design 
information  that  may  be  derived  from  the  numerical  simulation  of  back- 
mixed  combustion  processes. 

DISCUSSION 

leothamal  Flm> 

The  results  of  the  momentum  and  mass  transport  studies  identify 
several  important  characteristics  of  the  turbulence  subaodcls.  The 
turbulent  viscosity  from  the  two-equation  submodel  varies  through  the  flow- 
field  as  is  physically  expected.  The  two-equation  suoaodel  demonstrates 
the  generation  and  dissipation  of  turbulence  kinetic  energy  near  the 
Jet  exit,  the  stagnation  point,  and  near  the  boundaries.  In  contrast,  the 


28 


r 


viscosity  from  the  aJr.ehralc  subnodel  is  spatially  unifona  throughout 
tlic  flowfiold  for  lsotlieraal  flow,  and  varies  only  with  a change  in  the 
kinc-tic  energy  of  the  inlet  streams.  The  impact  of  these  differences 
between  the  two  eddy  viscosity  submodels  is  evidenced  by  the  velocity 
profiles,  the  predicted  location  of  the  stagnation  point,  the  radial  mass 
transport,  and  the  dissipation  of  the  kinetic  energy  of  turbulence  in  the 
wake  region. 

Velocity  Profiles.  The  predicted  velocity  profiles  agree  qualitatively 
with  experimental  results.  The  two-equation  submodel  more  closely 
approximates  the  experimental  trends  although  the  predicted  velocity 
contours  for  conditions  of  strong  recirculation  are  quantitatively 
incorrect  throughout  the  flowfield  for  both  turbulence  submodels.  The 
departure  from  experiment  is  especially  evident  near  vails,  and  suggests 
the  need  to  refine  the  wall  functions  used  to  specify  boundary  conditions 
for  turbulent  kinetic  energy  or  dissipation  rate,  and/or  to  decrease 
the  grid  spacing  near  the  wall. 

Stagnation  Point.  The  location  of  the  stagnation  point  as  indicated 
by  the  velocity  profiles  and  upstream  extent  of  tracer  mea'suroments  is 
better  described  by  the  two-equation  submodel.  The  effective  viscosity 
predicted  by  the  algebraic  submodel  is  spatially  uniform  throughout  the  flow- 
field  for  given  inlet  conditions.  No  account  is  cade  for  changes  in  the 
kinetic  energy  of  turbulence  or  energy  dissipation  rate  in  the  region  of 
encounter  between  the  main  and  jet  flows.  The  lower  viscosity  predicted 
in  tha  recirculation  zone  by  the  two-equation  submodel  is  more  effective 
in  describing  the  actual  axial  momentum  and  taass  exchange. 

Radial  Mass  Transport.  Although  the  experimentally  determined  tracer 
isopleths  arc  in  closest  agreement  in  the  recirculation  zone  for  the  two- 
equation  submodel,  the  downstream  radial  spread  in  more  effectively 
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described  by  the  algebraic  submodel.  The  uniformly  high  viscosity 
predicted  by  the  algebraic  turbulence  submodel  aids  the  overall  mass 
transport  throughout  the  flowfield,  i.e.  axial  and  radial  diffusion  from 
the  jet  is  greater  for  the  algebraic  rather  than  for  the  two-equation  sub- 
model . 

The  radial  mass  transfer  predicted  by  the  two-equation  turbulence 
submodel  suggests  that  the  nrlncioal  mechanism  of  mass  transport  is  hy 
large-scale  convection  rather  than  by  snail  scale  gradient  diffusion 
processes.1®  The  tracer  concentration  varies  gradually  through  the 
recirculation  zcr.e.  The  downstream  mixing  layer  is  bounded  by  steep 
concentration  gradients  that  enclose  the  jet  mixing  region.  The  numerical 
predictions  indicate  that  at  high  approach  velocities  the  bulk  flow  detours 
arour.d  the  recirculation  core,  retaining  upstream  flow  properties  (low 
viscosity  and  small-scale  mixing)  and  effectively  precludes  radially- 
directed  mass  transport.  Conversely,  the  experimental  trend  reinforces 
the  conjecture  that  turbulent  exchange  processes  should  be  active  • . 
downstream  from  the  recirculation  zone. 

Wake  Region.  The  kinetic  energy  of  turbulence  predicted  by  the 
two-equation  submodel  is  quickly  dissipated  in  the  wake  region  downstream 
of  the  recirculation  rone.  The  resulting  low  .viscosity  impedes  turbulent 
mixing  throughout  the  wake  region.  In  contrast,  the  uniformly  high 
viscosity  of  the  algebraic  submodel  is  effective  in  predicting  the 
radial  spread  of  tracer  observed  experimentally. 

Sot  Fla) 

The  hot  flow  results  reflect  the  inherent  limitations  of  the  under- 
lying eddy  viscosity  submodels  and  emphasize  the  urgency  to  refine  the 
basic  fluid  dynamic  representations  employed  in  the  analytical  models. 
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In  general,  the  predicted  hot  flow  properties  qualitatively  simulate 
gross  flame  features.  Before  adequate  experimental  correlation  is 
achieved,  however,  the  submodels  of  eddy  viscosity  and  turbulent  mass 
exchange  must  perform  satisfactorily  for  the  case  of  isothermal  flow, 
and  the  chemical  submodels  and  the  coupled  turSilence/chcnistry  interaction 
must  be  fully  explored  and  likely  improved  for  the  case  of  hot  flow. 
Refinement 

The  results  above  identify  directions  for  refinements  to  the 
eddy  viscosity  submodels,  inlet  conditions,  and/or  surface  boundary 
conditions.  A parametric  investigation  vas  initiated  to  inprove  the 
predictions  of  the  two-equation  eddy  viscosity  submodel  by  exploring 
the  effect  of  boundary  condition  specifications. 

The  outcome  of  the  parametric  studies  suggested  that  in  confined 
flows,  where  the  recirculation  zone  may  contact  solid  boundaries,  the 
specification  of  near  wall  turbulence  energy  dissipation  rate  (c)  may 
significantly  influence  the  nixing  characteristics  of  the  two-equation 
turbulence  submodel.  These  results  ore  consistent  with  the  findings  of 
earlier  studies11  where  a reduction  of  c adjacent  to  a downstream  facing 
wall  was  used  to  improve  numerical /experimental  correlation.  Similarly, 
the  reduction  of  € cn  the  upstream  facing  step  of  the  jet  tube  was  found 
to  radically  alter  the  predicted  mass  transport  behavior  of  the  two- 
equation  submodel.  The  Improved  radial  mass  transport  is  evidenced  by 
the  tracer  concentration  profiles  in  Figure  14  with  a slight  sacrifice  in 
stagnation  point  correlation.  The  stagnation  point  correlation  can  be 
improved  by  modifying  k or  c at  the  Jet  inlet. 
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The  refinements  clearly  demonstrate  the  sensitivity  of  the 
numerical  solution  to  boundary  condition  specification,  and  the  need 
to  test  the  two-equation  turbulence  submodel  under  isothermal  flow 
conditions  prior  to  hot  flow  simulations. 

CONCLUSIONS 

Numerical  predictions  of  the  turbulent,  backmixed  flow field  of  an 
opposed- Jet  combustor  have  been  compared  to  experimental  observations. 

Although  favorable  qualitative  correlation  has  been  established, 

| deficiencies  in  the  eddy  viscosity  submodels  and  associated  boundary 

| conditions  used  as  a foundation  for  the  transport  mechanisms  have  been 

| Identified. 

1 

* The  flowfield  hydrodynamics  are  more  accurately  represented  by  the 

| tvx>-cquation  submodel  because  of  the  intrinsic  variable  viscosity  and 

length  scale.  Comparable  solutions  may  be  obtained  at  considerably 
less  expense  by  employing  the  algebraic  submodel. 

» The  mass  transport  predictions  of  the  algebraic  submodel  demonstrate 

acceptable  correlation.  The  prediction  of  radial  mass  transport  by  the 
two-equation  submodel  in  general  does  not  confom  to  experiment.  The 
performance  of  the  tve-equation  submodel 'may  be  improved  by  careful 
specification  of  boundary  conditions.  In  particular,  the  near  wall 
turbulence  energy  dissipation  rate  was  shown  to  significantly  influence 

| 

the  mixing  characteristics  of  the  twe-equation  submodel. 

The  results  of  the  isothermal  flow  studies  indicate  that  the 
Tilgcbraic  turbulence  submodel  performs  adequately  in  describing  general 
fluid  flow  patterns  and  mass  transport  for  the  conditions  investigated. 
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A major  impact  of  the  present  study  ir.  that  considerable  testing  of 
the  numerical  model  is  required  for  isothermal  flow  before  processing 
to  the  complicating  conditions  of  combustion.  The  testing  requires 
a coupled  numerical/cxpcrimcntal  study  and  is  likely  governed  by  the 
geometry,  inlet  and  solid  boundary  conditions  of  the  system  under 
investigation.  The  hot  flow  calculations  presented  here  accentuate 
the  inadequate  transport  characteristics  identified  in  the  isothermal 
flow  analysis.  Attempts  to  refine  the  chemistry  submodel  or  to  quantify 
the  Interaction  of  fluid  motions  on  chemical  reaction  rates  are 
dependent  on  the  correct  turbulence  model  formulation. 


3. 4 HO  Probe 


Effects  of  oxygen  containing,  and  carbon  monoxide  contain- 
ing atmospheres  on  nitrogen  oxides  concentration  were  evaluated 
during  the  current  reporting  period.  Sample  probe  materials 
tested  included  304,  316,  and  321  stainless  steel  and  silica* 

The  length  of  each  sample  probe  was  arbitrarily  chosen_ 
to  be  2 meters.  The  residence  tine  of  the  doped  carrier  gas  in 
the  sample  probe  test  section  was  approximately  1 second  for 
the  4 litcrs/minute  flow  rate  and  2 meter  sample  probe  length. 

The  NO  and  NO^  input  levels  to  the  sample  probe  (Sample 
Point  2 in  Figure  7)  were  chosen  to  be  SOOppm  and  75ppn  res- 
pectively. These  levels  simulate  NO  concentrations  which  are 
typically  encountered  in  sampling  and  a 1I02  concentration  which 
is  thought  to  be  representative,  CO  levels  selected  were  0, 

100,  1000,  and  2500ppm. 

The  results  arc  presented  in  Figure  15  for  mixtures 
without  carbon  monoxide  and  in  Figure  16  for  mixtures  with 
carbon  monoxide.  The  percent  change  of  NO  and  N02  represent  the 
percent  change  in  concentration  between  Sample  Points  2 and 
3 (Figure  7^ . Since  total  NO^  was  conserved  in  all  cases,  the 
figures  demonstrate  a direct  correspondence  between  NO  and  N02 
transformations.  The  temperature  shown  is  the  gas  temperature 
at  Sample  Points  2 and  3.  The  results  without  carbon  monoxide 
were  presented  at  the  First  International  Chemical  Congress  in 
the  Western  Hemisphere  in  December,  1975. 52  The  results  with 
carbon  monoxide  were  presented  at  the  Western  States  Section  of 


the  Combustion  Institute  in  April,  1576. ls. 
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RESULTS  (WITHOUT  CARBON  MONOXIDE)  ' 
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The  results  fcr  silica  arc  presented  in  Figure  15a.  No  1 
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significant  transf ornation  is  observed  over  the  tenperature  \ 4 

m 

range  and  the  residence  tiae  studied.  Slight  oxidation  of  j j 

1 

NO  to  NO2,  proportional  to  oxygen  content,  is  discernible.  j t 

n 

The  304  stainless  steel  probe  (Figure  15b)  catalytinally  ► 

I 

reduces  HOj  to  NO  at  temperatures  in  excess  of  100*C.  The  | 

MSB 

conversion  of  NO^  to  HO  at  elevated  tenpcraturcs  is  consistent  i 

1 

with  the  results  of  a variety  of  studies  in  which  stainless  £ 

•1 

steel  has  been  evaluated  for  use  as  the  converter  waterial  in 

B 

choBiluninesceat  nitrogen  oxides  analyzers.1*  11  At  tenpertureo  : J 

m 

below  the  catalytically  active  tenperature  of  200“C,  no  signifi-  ' 9 

1 

cant  change  is  observed.  The  oxygen  content  of  the  mixture  has  3 

i 

an  important  impact  on  the  conversion  efficiency  at  the  peak  ‘*9 

E 

Veopcrature  evaluated  (400cC).  J9 

B 

The  results  fcr  *'he  .316  stainless  -steel,  presented  in 

m 

Figure  15c,  indicate  tha.  a higher  temperature  (>300*C)  is  ^9 

9 

r«  quired  to  initiate  the  catalytic  reduction  of  N<>2  in  coapar- 
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ixon  to  304  stainless  steel.  At  400'C,  however,  dissociation 

jB 

is  noar’y  1002  complete.  The  oxygen  content  of  the  mixture  has  (9 

i 

a nodest  impact  on  the  conversion  efficiency  at  the  peak  tenp-  41 

E 

crature  evaluated  (400*C).  Slight  evidence  of  oxidation  of  HO  29 

to  KOj  in  again  discernible  at  the  lower  teaperatures.  99 
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The  results  for  the  321  stainless  steel,  presented  In 
Figure  15d,  show  behavior  sinilar  to  the  316  stainless  steel. 

The  three  stainless  steels  evaluated  have  similar  carbon, 
chromium,  and  nickel  con  ten t . 1 7 ~l  9 As  shown  in  Table  VI, 
however,  316  stainless  steel  has  molybdenum  added  to  improve 
corrosion  resistance,  and  321  stainless  steel  has  titanium  add- 
ed to  reduce  carbide  precipitation.  In  general,  stainless  steels 
have  strong  resistance  to  corrosion  due  to  an  oxidized  layer 
on  the  surface.  Reducing  atmospheres  can  remove  the  oxide 
film  and  change  the  passivity.  It  is  likely  that  prior 
exposure  of  stainless  steel  probes  to  varying  sample  environ- 
ments may  alter  the  degree  of  NO^  transformations.  Evidence  of 
probe  history  effects  have  been  reported  in  the  literature,26 
but  a full  assessment  has  yet  to  be  made. 


RESULTS  (WITH  CARBON  MONOXIDE) 


The  results  for  silica  are  presented  in  Figure  16a.  Ho 
significant  transformation  is  observed  at  temperatures  below 
300*C.  The  slight  oxidation  of  NO  to  NOj  observed  without 
carbon  monoxide  is  not  discernible  when  carbon  monoxide  is 
present.  At  400®C  some  reduction  of  HO2  to  MO  occurs.  The 
reduction  is  c ndent  on  both  oxygen  content  as  well  as  carbon 
monoxide  At  the  level  of  carbon  monoxide  for  which  the 

transformation  is  most  pronounced  (2500ppm  CO),  the  dependency 
"on  ox>gen  is  the  reverse  of  the  dependency  at  lower  carbon 
monoxide  levels  (100,  lOOOppc  CO). 
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The  results  for  316  stainless  steel  are  presented  In 
Figure  16b  (304  and  321  stainless  steels  were  not  evaluated). 

The  presence  of  carbon  nonoxide  effectively  lowers  the  tern-* 
perature  at  which  reduction  of  N02  to  HO  is  observed.  Signif- 
icant reduction  is  observed  at  tenperatures  in  excess  of  100#C. 
The  reduction  is  more  pronounced  as  carbon  monoxide  is  increased. 
Increased  oxygen  content  has  the  effect  of  reducing  the  conver- 
sion at  a given  temperature  and  a given  carbon  monoxide  level. 


DISCUSSION 


The  results  without  carbon  monoxide  are  expected  and 
explainable  by  the  known  catalytic  performance  of  stainless 
steel  in  reducing  N02  to  The  results  with  carbon 

monoxide  are  new  and  yill  require  further  testing  to  fully 
assess  the  controlling  mechanism.  Nevertheless,  tentative 
explanations  are  possible  by  considering  the  competition 
between  the  following  two  reactions: 


C0  + 2°2  ' C02  • (S) 

catalyzed 

N02  + CO  HO  + C02  (10) 

Increased  CO  will  advance  the  reduction  of  K02'  by  reaction  (10) 
as  observed  experimentally.  Increased  02  will  Inhibit  the 
advance  of  reaction  (10)  due  to  the  competition  between  Oj  and 
KO  for  the  CO.  In  studico  of  CO  oxidation  by  HO  and  02,  oxygen 
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is  observed  to  preferentially  oxidize  CO  to  COj.’*'21 

The  silica  results  with  carbon  monoxide  indicate  that 
reaction  (10)  nay  be  active  at  iOO'C.  The  advance  of  reaction 
(10)  as  a homogeneous  reaction  at  400°C  is  slow  relative  to 
the  one  second  residence  tine.  The  dependency  upon  oxygen 
is  not  explained. 

The  results  of  Figure  16b  (316  stainless  steel  with  carbon 
nonoxide)  differ  from  an  earlier  investigation  by  Halstead2' 
in  t‘iat  total  oxides  of  nitrogen  are  conserved  .in  the  present 
study.  The  removal  of  NO  ovserved  by  Halstead  nay  result  from  5 

(1)  the  reduction  of  NO^  by  reducing  species  other  that  CO  j 

(e.g.  hydrocarbons)  present  in  the  fuel  rich  flames  sampled,  f 

(2)  the  elevated  (S  800*C)  temperatures  to  which  the  samples  i 

were  temporarily  exposed  at  the  probe  inlet,  (3)  the  longer  I 

residence  tine  (4  seconds  compared  to  1 second)  or  (4)  a com-  5 

bination.  Further  study  that  addresses  additional  reducing  ‘ 

species,  higher  temperatures  (>  400”C),  and  longer  residence 

times  (>  1 second)  will  be  necessary  to  fully  assess  the  difference. 


CONCLUSIONS • 

• Nitrogen  oxides  are  sampled  from  combustion  sources 
for  a variety  of  reasons  and  over  a wide  range  of 
temperatures  and  gas  composition. 

• A number  of— reactions  arc  potentially  active  to 
transform  nitrogen  oxides  concentration  while  sampling 
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combustion  products. 

• HO^  is  conserved  in  silica  and  304/  316,  and  321 
stainless  steel  sample  probes  for  the  temperature 
range  25  to  400*C  and  a residence  time  of  1 second. 

• In  the  absence  of  carbon  monoxide,  reduction  of  NO2 
is  observed  in  316  and  321  stainless  steel  at  temper- 
atures exceeding  300°C.  In  304  stainless  steel,  N(>2 
reduction  is  observed  at  temperatures  exceeding  100*C. 

• In  the  presence  of  carbon  monoxide,  the_  reduction  of 
HO2  occurs  at  temperatures  exceeding  100*C  in  316 
stainless  steel, 

• In  silica,  no  significant  transformation  of  nitrogen 
oxides  is  observed  in  the  absence  of  carbon  nonoxide. 

When  carbon  monoxide  is  present  'at  levels  exceeding 
lOOOppm,  reduction  of  NO^  is  observed  at  temperatures 
exceeding  300*C. 

• Previous  KO^,  NO  and  NO2  data  collected  at  moderate 
temperatures  (25  to  600*C)  with  stainless  steel 
sampling  lines  and  at  temperatures  above  200*C  with 
silica  must  be  used  with  caution. 

• Tbe  potential  for  chemical  transf ormatlon  cf  nitrogen 
oxides  concentration  while  sampling  combustion  products 
at  moderate  temperatures  may  be  reduced  by  (1)  using  only 
silica,  (2)  maintaining  the  temperature  below  300*C  in 
carbon  monoxide  containing  atmospheres,  and  (3)  exercis- 
ing appropriate  precaution  In  handling  water  lu  the 
sample 
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• Additional  information  Is  needed  to  assess  the  extent 
of  chemical  transformation  that  occurs  under  the  variety 
of  sampling  conditions  experienced  In  practice,  and* 
guide  the  selection  of  probe  and  sample  line  materials 
in  order  to  prevent  the  occurrence  of  chemical  transfor- 
mation. 


3.5  Publlcations/Papcrs  Presented 

Two  publications  were  submitted  during  the  reporting  period 
and  three  papers  were  presented.  Papers  are  now  in  preparation 
or  planned  for  publication  or  presentation.  The  current 
status  of  the  publications  and  presentations  of  papers  is 
summarized  in  Table  VIZ. 
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4.  DESCRIPTION  OF  RESEARCH  DIRECTION;  * 


CONTINUATION  YEAR 


I 


4.1  Model  Development  and  Evaluation 

In  the  continuation  year,  the  emphasis  will  be  placed 
on  the  evaluation  of  the  predictive  methods  for  the  case  of 
heated  (nonreacting)  flow  and  hot  (reacting)  flow.  The  sensi- 
tivity study  for  initial  and  boundary  conditions  for  the  case 
of  isothermal  (nonreacting)  flow  will  be  completed  as  well. 

Heated  Flow,  The  heated  case  is  designed  to  test  the  coupled 
submodels  of  eddy  viscosity,  mass  transport,  and  energy  transport. 
The  heated  flow  case  was  initiated  during  the  current  reporting 
period  with  the  design  and  installation  ,of  a preheat  system  into 
the  combustion  facility. 

Hot  Flow.  The  hot  flow  case  will  substantially  increase  the 
complexity  of  the  evaluation.  The  added  considerations  include 
chemistry  submodels  (including  the  effects  of  turbulence  on 
chemical  reaction  rates),  and  radiation  submodels.  Predicted 
profiles  of  velocity,  turbulence  kinetic  energy,  equivalence 
ratio,  temperature  and  concentrations  of  hydrocarbons,  oxygen, 
carbon  monoxide,  nitric  oxide,  nitrogen  dioxide,  snd  carbon 
dioxide  will  be  compared  to  experimental  results.  Comparison 
of  the  predicted  distributions  to  those  determined  experimentally 
will  be  used  to  identify  deficiencies  in  the  models  of  turbulence 
and  kinetics  for  the  condition  of  reacting  flow. 

The  hot  flow  case  was  initiated  during  the  past  year  by 
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(1)  generating  first  solutions  for  the  hot  flow  case  (Figures  12 
and  13),  (2)  generating  hot  flow  experimental  data  for  mcthanc/air , 
and  (3)  preparing  and  submitting  a proposal  to  the  National 
Science  Foundation  for  support  to  complement  the  hot  flow 
evaluation  studies. 

4.2  Mechanisms  of  Pollutant  /Plume  Formation 

The  experimental  studies  designed  to  specifically  explore 
the  mechanism  of  pollutant  and  plume  formation  will  commcnse 
in  the  summer  of  1976.  The  experiment  will  be  used. alone  and 
In  coabinaticn  with  numerical  results  in  the  next  twelve  month 
period  to  explore  pollutants  and  plume  formation.  Parametric 
studies  will  be  conducted  to  establish  the  relative  effect  of 
the  transport  ar*  chemical  processes  (chemical  kinetics,  mass 
and  heat  diffusion,  and  fluid  mechanics)*  and  system  parameters 
(jet  and  main  stream  flow  rates,  composition,  temperature,  size 
and  shape  of  the  recirculation  zone)  on  pollutant  formation. 

4.3  Supplemental  Studies 

Under  supplemental  studies,  the  exploration  of  the  lean 
stabilized  combustion  performance  of  the  OJC  will  be  initiated  and 
completed,  the  air  prehea.t  studies\will*  continue,  a nonpremixed 
combustor  configuration  will  be  designed  and  tested  to  complement 
the  prcnlxed  OJC  geometry,  swirl  will  be  added,  the  NO  probe 
study  will  be  expanded  to  explore  a variety  of  hydrogen  containing 
species  and  probe  history,  and  studies  will  be  initiated  that 
address  particulate  formation  and  luminosity,  especially  in 
relation  to  the  on-going  alternative  fuel  studies  ct  AFAPL.  The 


supplemental  studies  have  been  added  to  specifically  address 


current  needs  of  the  Air  Force. 
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5.  SUMMARY 


I 

j 


ii 

The  UCX  investigation  of  exhaust  pollutant  and  plume  forma- 
tion in  continuous  combustion  is  a combined  analytical  and 
experimental  study  of  turbulent,  backmixed  combustion.  Predicted 
profiles  of  the  flow  properties  are  being  systematically  compapcd 
to  experimentally  determined  profiles.  Model  development  and 
evaluation,  and  studies  of  the  mechanisms  of  pollutant  and  plume 
formation  are  being  conducted  for  premixed  methane/air  and  propane/ 
air  reactants. 

Future  work  will  include  nonpremixed  combustion  and  a 
liquid  fuel  system. 


5.1  Model  Development  and  Evaluation 

Nonroacting.  The  isothermal  evaluation  of  the  PISTEP  II 

; and  TEACH  numerical  codes  has  addressed  the  effective  viscosity 

■ and  mass  transport  submodels.  Still  to  be  completed  is  a thorough 

| investigation  of  initial  and  boundary  conditions. 

"he  heated  flow  evaluation  has  been  initiated  by  the  design  «• 

and  Installation  of  preheat  capability.  The  evaluation  will 

conclude  during  the  continuation  year. 

In  cooperation  with  the  EFA,  the  CRISTY  numerical  code  is 

being  run  as  tine  permits  in  parallel  with  the  PISTEP  II  and 

TEACH  numerical  codes.  The  interaction  is  proving  useful  to  both 

parties  (UCI/AFOSR  and  EPA),  and  comparati  e results  will  likely 

be  "available  ■ for  presentation  in  the  continuation  year. 

* 

Reacting.  Hot  flow  solutions  have  been  obtained  with  TEACH 
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and  the  numerical  codes  (TEACH  and  PISTEP  II)  are  currently 
being  expanded  to  include  propane  kinetics.'  Initial  experimental 
data  have  been  obtained  in  preparation  for  the  coupled  numerical/ 
experimental  evaluation  in  the  continuation  year. 

5.2  Mechanisms  of  Pollutant /Plume  Formation 

The  experimental  facility  is  undergoing  final  modification 

(inclusion  of  preheat  capability)  prior  to  initiating  the  scheduled 
exploration  of  the  mechanisms  responsible  for  pollutant  and 
plume  formation  in  continuous  combustion.  To  fjcproye  the  signal 
aquisition  of  the  LDV  System,  the  addition  of  a counter  mode 
to  complement  the  existing  frequency  tracking  mode  is  under 
investigation.  The  sampling  performance  of  the  combustion  test 
facility  is  being  improved  as  a result  of  the  supplementary, 

HO  probe  studies. 

5.3  Supplementary  Studies 

The  NO  probe  study  has  shown  significant  reduction  of 
nitrogen  dioxide  (NO2)  to  nitric  oxide  (NO)  can  occur  in  sample 
probes  and  sample  lines,  especially  in  the  presence  of  carbon 
monoxide  (CO)  at  modest  temperatures  in  stainless  steel.  The 
effects  of  probe  history  .and  reducing  species  in  addition  to  CO 
will  be  evaluated  in  the  continuation  year.  In  addition,  lean 
combustion  emission,  air  preheat,  non-prenixed , swirl,  and 
particulate/luminosity  studies  will  be  initiated. 
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TABLE  I 

Elliptic  Flow  Codes 

CODE 

PRINCIPAL 

VARIABLES 

OR1CIKATOR 

REFERENCE  - • 

PXSTEP1 1 

^ - u 

laperial  College 
London 

3,25 

TEACH 

V - v 

laperial  College 
London 

26 

CR1STY* 

* - » 

Dotted  Technology 
Research  Center 
Connecticut 

27,28 

* Earlier  vcrloas  called  FKEP  and  CFREP 
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TABLE  IV 


Predictive  Methods 


Xaproveaent s/Mod if lest  loos  During  Reporting  Period 


following  is  a description  of  tasks  which  required  either 
•odlf .'cation  of  existing  methods  or  new  developments 


TFACH 

1.  lsotherosl  flow  - submodels  for  notstntua  and  case 
transport  totted  against  experiment.  Itsproved  boundary 
conditions  required. 

2.  Hot  flow  - solution  obtained  for  methane  oxidation 
including  foraation  of  nitric  oxide. 

PISTEP  II 

1.  lsotherosl  flow  - tests  oj  submodels  for  acncntun  and 
maos  transport  against  experiment  initialed.  Evidence 
of  boundary  condition  modification  observed. 

CRISTT 

1.  Isotheraal  flow  - successfully  modified  for  VCJ  OJC 
geometry  in  cooperation  with,  the  EPA.  Tests  against 
experiment  initiated. 


Following  1*  * description  of  remaining  problems  and  plans 
for  their  resolution. 


TEACH.  PISTEP  11.  CBtSTY 

Mo  problems  are  current  excluding  the  chronic  demand 
of  numerical  work  for  additional  resourcee.  Ths 
development  and  evaluation  of  the  numerical  codes  arc 
modestly  behind  schedule  t but  the  present  rete  of  progress 
and  success  in  securing  additional  resources  will  return 
the  development  and  evaluation  to  the  planned  schedule. 


TABLE  V 


Combustor  Test  Facility 

leprovenen t s/Kodif i c«* t ion*  During  Kcportlng  Period 


Following  is  a description  of  tasks  which  required  cither 
modification  of  existing  equipment  or  new  work  to  be  performed. 

1.  Inlet  air  flot?  preheat  - to  core  effectively  simulate  gas 
turbine  combustion  and  explore  the  effect  of  air  preheat 

on  nuuurieal  model  development  and  evaluation  and  pollutant/ 
plurte  formation,  air  preheat  capability  vac  added.  Inlet 
temperatures  to  600*C  (1122CF)  are  nov  available. 

2.  Exhaust  heat  load  - an  air  dilution  system  and  vater  cooled 
heat  exchanger  coil  vcrc  installed  to  allov  higher  combustor 
loadings  to  be  accomodated  without  thermally  taxing  the 
building  exhaust  system. 

3.  Sample  probe  - a 1 /8-inch  O.D.  sample  probe  vat  designed , 
built,  tested , and  used  for  composition  measurements. 

The  1/8-inch  size  reduced  the  impact  on  the  flow field  of 
the  previous,  1/4-inch  probe. 

4.  Emission  Console  overhaul  - the  emission  console  was 
thoroughly  inspected  and  serviced  prior  to  the  generation 
of  the  data  base  rune.  All  tubing  was  disconnected, 
cleaned,  and  reassembled. 

5.  Data  Acquisition  System  - the  DEC  PDP-11  data  acquisition 
system  was  installed  and  successfully  interfaced  with  the 
combustor  facility  instrumentation.  All  necessary  software 
has  been  developed,  cheeked  and  is  currently  in  operation. 


Following  is  a description  of  remaining  problems  and  plans 
for  their  resolution. 

1.  Combustor  Site  - the  Sinn  (2-inch)  J.D.  combustor  site 
has  proven  edequate  for  the  studies  conducted  but  not 
optimum  for  the  resolution  desired  for  property  profile 
measurements.  A 78a t (3-inch)  combustor  section  has  been 
design'd  and  is  currently  being  fabricated ,\  The  flexibility 
of  changing  combustor  size  will  prove  useful  in  evaluating 
boundary  condition  inpael  on  the  numerical  codes. 

2.  Sample  probe  design  - etudiee  conducted  to  explore  chemical 
transformation  of  nitrogen  oxides  in  probes  and  lines 

have  shown  that  stainless  steel  will  eatalytieally  reduce 
*°jk  to  SO  in  the  pretsr.ee  of  CO  at  temperatures  in  excess  _ 
of  )00*C.  A silzes,  water  coated  sampling  probe  is 
equipped  with  aerodynamic  quench,  is  currently  being 
designed. . 

3.  Swirl  - awirl  ie  an  integral  factor  in  promoting  plume 
stabilization  in  practical,  continuous  combustion  devices. 
The  numerical  methods  must  be  tested  for  the  case  of  twirl 
es  well  as  non-twirl  if  successful  simulations  of  practical 
devices  arc  to  be  achieved.  Swirl  vanes  are  currently 
being  designed. 

4.  LOT  System  - the  perfcrmcr.ee  of  the  frequency  tracker 
deteriorates  lo  ucaccrpicMc  levels  in  the  highly 
turbulent  region  of  Jet  penetration.  Seeding  techniques 
for  the  Jet  and  expansion  of  the  signal  acquisition  system 
tc  include  the  counter  mode  arc  bring  explored. 
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